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Abstract 

, We propose an effective theory which governs Pomeron dynamics in QCD at high energy, in the 

leading logarithmic approximation, and in the limit where N c , the number of colors, is large. In 
spite of its remarkably simple structure, this effective theory generates precisely the evolution 
, equations for scattering amplitudes that have been recently deduced from a more complete 

microscopic analysis. It accounts for the BFKL evolution of the Pomerons together with their 
r S^ • interactions: dissociation (one Pomeron splitting into two) and recombination (two Pomerons 

merging into one). It is constructed by exploiting a duality principle relating the evolutions in 
the target and the projectile, more precisely, splitting and merging processes, or fluctuations 
in the dilute regime and saturation effects in the dense regime. The simplest Pomeron loop 
calculated with the effective theory is free of both ultraviolet or infrared singularities. 
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There has been recently significant progress in our understanding of high energy 
hadronic scattering, and in particular of the processes occurring at large parton densities 
and which are believed to be responsible for the unitarization of the scattering amplitudes 
and the saturation of the parton distributions. Non linear evolution equations have been 
derived which describe the approach towards saturation and the unitarity limit, and which 
have the structure of stochastic evolution equations. However, it has been very recently 
recognized [1] that the equations which were considered as the most complete, namely the 
Balitsky-JIMWLK (Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner) equa- 
tions [2-5], are in fact incomplete. This is manifest in the statistical language by the 
presence of fluctuations at high momenta [6, 7] which are not well accounted for by the 
JIMWLK evolution of the target wavefunction [1]. In the language of Pomerons, the 
JIMWLK equation contains Pomeron merging but not also Pomeron splitting 2 . 

Following this observation, two of us (E.I. and D.T.) have constructed a hierarchy of 
non-linear evolution equations for the dipole scattering amplitudes which include both 
gluon mergings and gluon splittings, and thus generate Pomeron loops through iterations 
[1, 8]. These equations have been argued to hold in the limit where the number of colors 
N c is large (N c ^> 1), and indeed it has been checked explicitly in Ref. [8] that the 
vertices appearing in these equations are the same as the corresponding 'triple Pomeron 
vertices' computed in perturbative QCD at large N c [9-11]. A complementary approach 
has been developed by Mueller, Shoshi, and Wong [12] who proposed a generalization of 
the JIMWLK equation which includes the effects of pomeron splitting in the dilute regime 
and for large iV c . These two approaches follow the same general strategy — namely, they 
combine the non-linear JIMWLK equation at high density with the color dipole picture 
[13, 14] in the dilute regime — and lead indeed to the same evolution equations for the 
scattering amplitudes, as demonstrated in Ref. [8]. (See also Refs. [15, 16] for related recent 
developments.) 

It is our purpose in this letter to show that the equations obtained in [1,8] and [12] can 
be reformulated in term of an effective theory for Pomerons. By 'Pomeron' we mean here 
the color singlet exchange which describes the interaction between an elementary color 
dipole and the field of a target in a single scattering approximation, and which reduces to 
two gluon exchanges in lowest order perturbation theory. The construction of the effective 
theory involves a projection onto restricted degrees of freedom, precisely the Pomerons, 
and is expressed in terms of a simple Hamiltonian which describes the BFKL evolution of 
the Pomerons together with their splitting and merging. By requiring that the evolution 
should lead to identical results whether it is viewed as the evolution of the target or that 
of the projectile, one arrives at a duality principle which is used to construct the effective 
Hamiltonian from the Hamiltonian derived in [12] in the dilute regime. The limitations of 
the effective theory, and the subtle mathematical problems that arises when one attempts 
to analyze its microscopic content will be briefly discussed at the end of this letter. 

2 Note that the opposite terminology for what one means by 'splitting' and 'merging' would be 
more natural in relation with Balitsky equations, which refer to the evolution of the projectile. 
To avoid confusion on this point, in this paper we shall systematically use the terminology 
appropriate to target evolution. 
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Most treatments of high energy scattering rely on an asymmetric approach: typically, 
the 'projectile' is viewed as a collection of test particles which probe the color field of 
the 'target'. At high energy, the eikonal approximation is a good approximation, and the 
scattering of an elementary color charge is described by a 'Wilson line' of the form 



V£[a] = Pexp (ig J dx~a a (x~, x)t a ^j , 



(1) 



where x denotes the transverse coordinate of the particle, which is not affected by its 
interactions with the field of the target a a (x~,x), t a are the generators of the SU(3) 
algebra in the representation appropriate for the test particle, and the symbol P indicates 
that, in the expansion of the exponential, the color matrices a a (x~ , x)t a must be ordered 
from right to left in increasing order in x~ (we are using light-cone vector notations, 
x ± = (t ± z)/y/2.). For a more complex projectile, viewed as a collection of elementary 
color charges, the S-matrix is given by a product of Wilson lines like Eq. (1), one for each 
elementary color charge. 

In a frame in which most of the total rapidity Y is carried by the target, the target 
wavefunction can be described as a color glass condensate [4,17], and the corresponding 
S'-matrix is obtained as: 



(S) Y = J D[a] W Y [a] S[a] , 



(2) 



where a = a a (x~ , x) is a classical field randomly distributed with weight function W[a] (a 
functional probability distribution), and S[a] is the projectile S'-matrix for a given con- 
figuration of this random field. With increasing Y, the weight function evolves according 
to a functional renormalization group equation, of the generic form: 



^W Y [a] = -H 



a, 



i8 a 



W Y {a], 



(3) 



where if is a functional differential operator commonly referred to as the 'Hamiltonian'. 
Alternatively, one can view the same evolution as a change in the scattering operator, 
for a fixed weight function W[a:]. To see that, take a derivative w.r.t. Y in Eq. (2), use 
Eq. (3), and perform an integration by part in the functional integral: 



-^-p (S) Y = ~ j DN W[a] H ] 



a 



' i5a 



S[a\. 



(4) 



This can be interpreted as describing the evolution of the scattering operator Sy[a], with 
'Hamiltonian' W: 



d_ 

dY 



S Y [a] = -H^ 



a. 



i5a 



S Y [a] 



(5) 
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Both points of view, somewhat reminiscent of, respectively, the Schrodinger and the 
Heisenberg pictures of quantum mechanics, will be used in the following discussion (al- 
though we shall refrain from introducing explicitly rapidity dependent operators). In the 
Schrodinger picture, one puts emphasis on the evolution of the state vector, whose role is 
played here by the weight functional WV[a]. In the Heisenberg picture, the state vector is 
a constant reference vector involved in the calculation of all expectation values, here W[a], 
and one puts all the evolution in the operators, here the scattering operators SV[ai]. The 
Schrodinger picture corresponds to evolution equations which aim at providing a detailed 
microscopic description of the color field in the target, together with its complicated corre- 
lations. This is what the JIMWLK equation does. The Heisenberg picture rather describes 
how the test particles get dressed by color field fluctuations as they are boosted to higher 
rapidities. In this approach, the complicated color correlations in the target wavefunction 
are not immediately visible, and indeed the resulting equation of motion are established 
somewhat more easily. This second approach is essentially the one used by Balitsky to 
obtain his hierarchy of equations. 

The test particles that we shall consider are in fact elementary color dipoles, whose 
scattering amplitude reads: 



T(x,y) = l-±tr(VlV v ), (6) 

for a dipole with the quark leg at x and the antiquark leg at y. Here the Wilson lines are 
taken in the fundamental representation. We shall be interested in situations where the 
dipoles scatter off the color glass in the two-gluon exchange approximation (weak field 
limit) and we shall work in a large-iV c limit. In the weak field limit, the amplitude for 
a single dipole to scatter is obtained after expanding each of the Wilson lines to second 
order in a: 

V%[a] = l + ig J dx~a a (x~,x)t a 

2 

- y J dx~J dy-a a (x-, x)a\y-, x) [9(x~ - y~)t a t b + 0(y~ - x-)t b t a 

+ (7) 

Note that, to this order, the x~-ordering of the color matrices starts to play a role in 
Eq. (7). Still, this ordering is irrelevant for the computation of the dipole amplitude to 
lowest order, because of the symmetry of the color trace: tr(t a t b ) = \5 ab = tv(t b t a ). 
Namely, one finds: 

T(x, y) ~ T (x, y) = ^ [a a (x) - a a (y)] 2 , (8) 

which involves only the integrated field a a (x) = J dx~a a (x~ , x). Similarly the amplitude 
for k dipoles to scatter is given, within the same approximation, by Tq K \xi, y\, x K , y K ) = 
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T (xi, y 1 )...T (x K , y K ). In what follows, we shall refer to the amplitude (8) describing the 
single-scattering of an elementary dipole off a given color field as to a "Pomeron ex- 
change". Similarly, describes the exchange of k Pomerons. 

At this point we find it useful to digress on the linear evolution equation known as 
BFKL equation. This will allow a few observations which illuminate some of the math- 
ematical subtleties involved in taking the large N c limit when constructing our effective 
theory. Consider first the JIMWLK Hamiltonian. As shown in Ref. [18], when it is re- 
stricted to act on gauge-invariant observables, it can be given the simple form: 



-^jimwlk — 

"16^ / M (^v,z)(l + VX-VlV z -V^V v ) 



ab 



u,v,z 



X iSa^(u) iSa b Y (vY (9) 



where Ai is the dipole kernel 



M(x,y,z)= {X V )2 (10) 
(x - zf{z - yf 

Here the Wilson lines are in the adjoint representation. The derivatives can be freely 
moved across the bilinear form in Wilson lines, because they commute with the latter 
in the presence of the dipole kernel. That is, -£Tumwlk is Hermitian. The BFKL limit of 
-^jimwlk is obtained by expanding the Wilson lines to lowest non-trivial order in a. One 
gets: 

2 

^bfkl = - T |^ / M(u,v,z)[a a (u)-a a (z)][a h (v)-a b (z)] 



u,v,z 



X f^f^ , cf \ X d ( \ ' (11) 

It is not difficult to verify that -£/bfkl is again Hermitian. 

Let us now turn to the 'large-iV c limit'. This is obtained by (i) restricting the action of 
-f^BFKL to the dipole operators mentioned above and (ii) preserving only the dominant 
terms at large N c in the action of the Hamiltonian on these operators. When acting on 
the color fields inside a single factor To (i.e., on the same dipole), the two functional 
derivatives in Hbfkl yield a factor 5 cd , and then f ac ff b f c = —N c 5 ab produces the expected 
N c enhancement. On the other hand, the action on the color fields within two different 
factors T (i.e., upon two different dipoles) produces no such enhancement. Thus, at large 
N c , -f^BFKL can be equivalently replaced by an effective Hamiltonian in which the two 
functional derivatives are traced over color. This Hamiltonian, which we denote Hq for 
reason which will become clear shortly, is 



c u,v,z 
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where a s = a s N c /-K. Let us emphasize that, as obvious from the construction we have 
given, the two derivatives in Hq are to act on the same dipole. Note also that, as opposed 
to the original Hbfkl, Hq is not Hermitian: in fact, it is readily seen that its adjoint 
is ill-defined. This reflects the fact that the construction of Hq involves a projection on 
a specific set of degrees of freedom, and once this is done, one looses the possibility to 
integrate by part as in Eq. (4) in order to let Hq act on the weight functional W[a]. 
These special mathematical properties, restriction of the space on which the Hamiltonian 
is acting and loss of hermiticity, are general, and peculiar, mathematical features of the 
effective theory that we shall present. It is tempting to speculate that in doing the large 
N c limit we are renouncing to follow the evolution of some color correlations (precisely 
those which are suppressed at large N c ). The corresponding loss of information may be 
responsible for the simpler Markovian stochastic theory that we shall arrive at. 

We now return to the main stream of our discussion and establish a useful property. 
In Ref. [19], a symmetric description was obtained for the scattering between two color 
glasses in the regime where both systems are in the weak field regime). The final formula 
reads 



(S) Y = J DM W Y - y [a R \ J D[a L ] W y [a L ] e'/^^t'KW . (13) 

In this expression, p a L (x) = —V 2 x a a L (x) is the classical color charge density of the left- 
mover, and WV-j^or] and Wy[ai~\ are the weight functions for the right-moving and, 
respectively, left-moving color glass (note that the rapidity of the left mover is measured 
positively to the left, so that as we vary y, the total rapidity interval between projectile 
and target remains equal to Y). The precise conditions for the validity of Eq. (13) are 
detailed in Ref. [19]. Let us emphasize here a non-trivial aspect of this formula. Although 
it is essentially a weak field formula which assumes that the elementary dipoles interact 
only once, it contains the possibility that any number of dipoles of the projectile interact 
with an equivalent number of dipoles in the target. Thus Eq. (13) does account for multiple 
scattering, albeit in a restrictive way (each dipole interacting only once). These multiple 
scattering generate unitarity corrections if Y is large enough. At the same time, we require 
both color glasses to be unsaturated. This imposes some limited range of variation for y 
within which Eq. (13) is correct. 

Now, Lorentz invariance implies that (S Y ) may depend on the total rapidity interval 
Y, but, within the range of validity of Eq. (13), cannot depend upon the rapidity y used 
to separate the system into a 'projectile' and a 'target', or equivalently on the frame 
which we choose to describe the collision. This implies (see also Ref. [20] for a similar 
argument): 



= 



d(S) 



Y 



dy 



^W Y -y[a R ]^W y [a L ] + W Y - y [a R ] W y [a L ] 



(14) 
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The evolution of both weight functions are given by: 



^ W Y -y[a R \ = -gy W Y - y [a R ] 



H 



OCR, 



i5a R 



W Y - y [a R ], 



^-W y [a L ] = -H 



S 



OIL, 



W y [a L }. 



(15) 



We shall keep the evolution of the left-mover as shown in the above equation, but perform 
an integration by parts in the functional integral over a R in Eq. (14). Next, we note that 



Using this identity in Eq. (14) and performing a further integration by parts, now w.r.t. 
Oi (recall that p%{x) = —^ 2 x ct a L (x)) , one is left with a differential operator acting on 
Hy[ax,] = Wylpi] (with a slight abuse in the notation): 



OCR, 



i5a 



rt. 



iSo, 



-,PL 



W y [ PL ]. (17) 

For Eq. (14) to be satisfied, the contribution above should cancel against the term in 
Eq. (15) describing the evolution of the left-mover. This condition leads to the 'self- 
duality' condition: 



iSor. 



,PL 



Wy[ PL \. (18) 

The same relation holds obviously for the 'right' variables a R , p R . 

Going back to Eq. (16), one sees that what is involved in the duality operation 3 is a 
matching of splitting processes in the left movers, encoded by terms in the Hamiltonian 
of the form p 2 5 n /5p n , into merging process in the right movers, corresponding to terms 
of the form a n S 2 /Sa 2 . An example of such a matching is illustrated in Fig. 1. Splitting 
terms dominate in the dilute regime where they control the fluctuations, while merging 
terms become essential in the saturation regime where parton densities are large. This 
fluctuation-saturation duality is turned into a constraint on the evolution Hamiltonian of 
either the projectile or the target in Eq. (18). 

The self-duality constraint, which we expect to hold within the limited range of ener- 
gies in which the factorization (13) is valid 4 [19], will be used now to construct a simple 

3 To our knowledge, the duality between the roles of the operators p 2 5 n /Sp n and a n 5 2 /5a 2 has 
been first recognized by L. McLerran. 

The self-duality condition (15) has recently been claimed to hold in a much broader context 
[20]. 



H 



OiL, 



i5ai 



W y [a L ] = 



i5or 



• PL 
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Fig. 1. An illustration of the relation (16). p denotes the color charge of the left -mover (L) and 
a is the color field of the right-mover (R). The same diagram (any of the two diagrams in the 
middle) can be viewed either as a merging in R, or as a splitting in L. The first interpretation 
is natural when the diagram is produced by acting on the eikonal line with the Hamiltonian 
H\_^ 2 ~ p 2 5 4 /5p' i for L (cf. Eq. (19)). The second interpretation rather corresponds to the 
action of ~ a A 5 2 /5a 2 for R (cf. Eq. (24). 

Hamiltonian describing the Pomeron dynamics in the dilute regime, starting from the 
known, dominant, contribution containing only splitting processes that has been con- 
structed by Mueller, Shoshi and Wong [12]. Quite remarkably, and somewhat unexpect- 
edly, this Hamiltonian leads to equations of motion which reproduces the exact ones at 
large N c [1,8], that is the effective theory appears to be valid beyond the dilute limit 
where it is established. The Hamiltonian constructed in [12] reads 

Q 2 OC s f 

Hi^ 2 = -J^p ^ J M ( u ^ v ^ z )^(ui\u, z)g(v l \u, z)g(u 2 \z, v)g(v 2 \z,v) 

x c w c ? \/ 2 \/ 2 a c (u)a c (v). (19) 

5a a {u±) Sa -^) 5a b (u 2 ) 8a b (v 2 ) 

In Eq. (19), the integration goes over all the transverse coordinates u, v, z, u±, Vi, u 2 , 
v 2 . The function Q(ui\u, z) is, up to a factor gt a , the classical field created at ui by the 
elementary dipole (u, z), and reads 

1 («i — z) 2 , , 

Q ( Ul \u,z) = — lnf^ '-. 20 

An (ui - u) 2 

It is easy to understand (and was explicitly shown in [8]) that this Hamiltonian gener- 
ates Pomeron splittings. More precisely, the result of the operation of H\^ 2 on the two- 
Pomeron exchange amplitude T is proportional to T , and thus generates the following, 
fluctuation, term in the evolution equation for Tq 2 \xi, yi, x 2 , y 2 ): 

H ^ T o 2) ={^) 2 ^ J M(u,v,w)A (x 1 ,y 1 \u,w)A (x 2 ,y 2 \w,v)W 2 v ylT (u,v), 

u,v,w 

(21) 

where a 2 Ao is the amplitude for dipole-dipole scattering in the two-gluon exchange ap- 
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proximation and for large N c 



A (x,y\u,v) = - 



In 



[x — V) 



(y - uf 



x — u) 2 (y — v) 2 



(22) 



Clearly, this process corresponds to the splitting of one Pomeron into two. In general the 
Hamiltonian H±^ 2 can describe the transition n — > n + 1, in which case n — 1 of the 
Pomerons are simply "spectators" . Note that H^ 2 is non-Hermitian, which we interpret 
as reflecting again the large-iV c approximation implicitly involved in its derivation. 



To apply the duality transformation, it is convenient to reexpress H\^ 2 m terms of 
the sources p a (x) of the color field a a (x), by using p a (x) = — V 2 x a a (x). We then obtain 



9 ®s_ 
IQNl 2tt 



J M(u,v,z) 



X 



u,v,z 

5 



1 2 



5p a {u) 5p a (z) 



1 2 



Sp b (z) 5p b (v) 



p c (u)p c (v). 



(23) 



At this point we force the Hamiltonian to be self-dual. This is done by adding to H^^b /i$P, p] 
its dual H{^ 2 [a.,5/i8a\ = H 2 ^\ (this new notation will be justified shortly). The Hermi- 
tian conjugate of H 2 ^i reads 



9 a 



2->l 



16N? 



^ / M(u,v,z)[a a (u)-a a (z)] 2 [a b (z)-a b (v) 
In J 



5a c (u) 5a c (v) 



, (24) 



and the action of H 2 ^ ± on the dipole scattering amplitude is 



HLi T (x, y) = ^jM(x, y, z) ^ [a a (x) - a a {z)f[a\z) - a\y)f 



27 



J M(x,y,z)T^\x,z;z,y). 



(25) 



Thus H 2 _^i generates the non-linear term in the first Balitsky equation. Similarly, it is 
obvious to show that the operation on Tq K \xi, y±; x K , y K ), will generate correctly the 
non-linear terms of the K-the Balitsky equation in the large-iV c limit (this is trivial; only 
one amplitude is "active", and we need to take into account all the possible permutations). 
Therefore the Hamiltonian in Eq. (24) generates in an effective way Pomeron mergings 
(hence the notation if 2 ->i); one has a transition of the form n + 1 — > n where, again, n — 1 
of the Pomerons are spectators. 



Thus the total Hamiltonian of our Pomeron effective theory reads 

= Hl + HU 2 + H ] 2 ^. 



(26) 
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The Hamiltonian Hq, describing the BFKL evolution, plays here the role of the free 
Pomeron Hamiltonian. The other two pieces i?|->i an d H\^2 correspond respectively to 
Pomeron merging and splitting, and will naturally generate Pomeron loops in the course 
of the evolution. The minimal Pomeron loop, which is simply the one-loop correction to 
the scattering amplitude {T(x,y)) Y , can be isolated by the successive operation of these 
two parts of the Hamiltonian, namely PL = /j|^ 2 ij]_^ 1 T . The explicit results reads 

PL = -g) 2 ^) 2 J M(x,y,z)M(u,v,w) (27) 



u,v,z,w 



x Ao(x,z\u,w)Ao(z,y\w,v)Viyi (T (u,v)) 



Y ■ 



Note that this result is free of any (ultraviolet or infrared) divergences. For instance, the 
poles in the dipole kernel at z = x and y are harmless because of Ao(x, z = x\u,w) = 0. 

A simple physical picture of this result is obtained by assuming that this Pomeron 
loop has been generated after the first two steps in the evolution starting with a target 
which is itself an elementary dipole (x ,y ). Then, Eq. (27) simplifies to: 



PL = -2 ( ^ 

.2-7T 



a x 2 



a 4 s J M{x,y,z)M(x ,yo,w)Ao(x,z\x ,w)Ao{z,y\w,yo). (28) 



This result has a clear physical interpretation: Both original dipoles — in the projectile 
and the target — split into new dipoles, processes which are represented by the two dipole 
kernels times a 2 s . Then, the child dipoles from the two systems scatter with each other, 
by exchanging two pairs of gluons; this yields the two factors A times a 4 . Finally, note 
that this contribution is negative, as expected, leading to a decrease in the amplitude in 
the course of the evolution. 

As we have already emphasized, the Hamiltonian (26) reproduces the complete equa- 
tions of motion established in [1,8] and [12]. While this intriguing property deserves 
further investigation, some insight can be gained by analyzing how the merging processes 
in the effective Hamiltonian compare to those deduced from correct microscopic dynamics 
as described by JIMWLK. The action of -£/jimwlk, Eq. (9), on the full dipole scattering 
amplitude T(x,y), Eq. (6), is 

XX 2 

T(x,y) = ^-(S yv - S xv ) \S ux tr(t b t a V^V y ) - 5 uy tr(t a t b V^V y )] (29) 



5a a {u) 5a b {v) v N t 
Simple algebra then easily yields the first Balitsky equation : 

H JlMWLK T(x, y) = -±J M(x, y, z)[-T(x, y) + T(x, z) + T(z, y) - T(x, z)T(z, y)\. 

(30) 

Then, after expanding the dipole operator T in the weak-field limit, and keeping terms up 
to the quartic order with respect to gauge field a, one finds an evolution equation which 
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contains not only the BFKL dynamics, but also the lowest order mergings (four gluons 
merging into two). 



Consider now the action of the JIMWLK Hamiltonian on T (x,y). Since: 



S S cp* 
8a-{u) 5a b (v) T °( X,y) = 2N C ~ S v»W™ ~ <W' 



(31) 



we have: 




^jM(x, y, z)Tr (l + V^V y - V^V Z - V}V y ) . (32) 



z 



When expanding the Wilson lines in the r.h.s. in powers of a, one obtains quadratic terms 
describing the BFKL evolution of T (x, y) plus higher order terms which describe n — > 2 
gluon mergings. But at this level, it is easy to see that the 4^2 terms generated by 
this expansion are not the same as those in the r.h.s. of Eq. (25). For instance, while the 
merging term in Eq. (25) includes a piece containing three different transverse positions 
(i.e., a^,a z a b z a b y ) , the corresponding JIMWLK result in Eq. (32) cannot generate such 



We thus see that thet actual, microscopic, dynamics of gluon merging in QCD is con- 
siderably more complicated than in our simple effective theory, yet the latter provides, as 
we have seen, the correct evolution equations for the scattering amplitudes. This shows 
that the additional merging terms generated by the JIMWLK Hamiltonian must compen- 
sate in the evolution equations against non-linear (quartic in a, or higher) contributions 
to the scattering amplitudes, as obtained by expanding the Wilson lines in equations like 
(6). We can refer to the latter as describing the dressing of the pomeron with multiple 
scattering. 

This brings us to comment on the nature of the dynamics described by the effective 
theory. This theory generates evolution equations for the Pomeron operators which 
are formally identical to the equations satisfied by the complete dipole scattering operators 
in QCD at large N c . This means in particular that the solutions to the equations 
for (T (K) ) r will appear to saturate the unitarity (or 'black disk') limit To = 1 in the high 
energy limit, in spite of the fact that the respective operators describe single scatterings 
only ! This indicates that one must be extremely careful in the physical interpretation of 
the effective theory. 

Let us then have a closer look at the microscopic dynamics that is describes. Effec- 
tively, the evolution of the target reduces to that of a system of dipoles subjected to a 
dynamics of a reaction-diffusion type: the dipoles undergo BFKL dynamics, they can split 
(one dipole into two dipoles), and they can also recombine with each other (two dipoles 
into one). The dynamics of such a system of dipoles is entirely coded in the /c-body den- 
sities riy' (see Sect. 5 in Ref. [1] for a precise definition). Although we shall not work this 
out explicitly here, it is not hard, by using the results of Ref. [19] to relate these dipole 
densities to colorless correlation functions of the color charge density p a . For instance the 




terms. 
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dipole number operator nix, y) can be identified with the bilocal operator p a (x)p a (y) of 
the effective theory. With such identifications, and by using Eq. (26), it is straightforward 
to construct the evolution equations satisfied by the dipole densities. One thus finds that 
ny(x, y) obeys the BFKL equation supplemented by a negative term proportional to ny\ 
which is generated by the merging piece i?2->i °f the Hamiltonian. Furthermore, the r.h.s. 
of the equation for drxy jdY includes the standard BFKL terms describing the individual 
evolutions of the two dipoles {xi,yf) and (2:2,2/2), but also a positive, fluctuation term, 
proportional to ny — this is generated by the splitting piece H\^ 2 of the Hamiltonian, 
and is the same as the corresponding term deduced from the dipole picture in Refs. [1,8] 

(3) 

- and, finally, a negative, recombination, term proportional to n Y . We thus obtain an 
infinite hierarchy, which describes a dipole reaction-diffusion dynamics, as anticipated, 
and predicts the saturation of the dipole density at a value of order l/a 2 s . 

Now, it is clear that this is only an effective dynamics since, as well known, dipoles in 
real QCD do not simply recombine with each other: the interaction between two dipoles 
inside the target wavefunction goes beyond the large-7V c approximation and leads to 
more complicated color configurations, involving higher color multipoles [13, 14]. The rea- 
son why it has been possible to simulate the non-linear effects responsible for unitarity 
corrections in the equations for the scattering amplitudes through simple 'dipole recom- 
bination' processes in the target wavefunction is because the same non-linear effects can 
be interpreed as projectile evolution, in which case they describe the splitting of a dipole 
in the projectile. Then, the 1 — > 2 dipole splitting vertex from the projectile is simply 
reinterpreted, within the effective theory, as a 2 — > 1 'dipole merging' vertex in the tar- 
get. Note finally that a similar dipole model including splitting and recombination has 
been recently used in Ref. [15] to generate evolution equations with Pomeron loops. The 
present work shows how this effective dynamics may indeed emerge from the actual target 
dynamics in QCD, and points to numerous subtleties involved in this precise connection. 
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